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ABSTRACT 

In  the  emerging  paradigm  of  open  spectrum  access,  cog¬ 
nitive  radios  dynamically  sense  the  radio-spectrum  environ¬ 
ment  and  must  rapidly  tune  their  transmitter  parameters  to  ef¬ 
ficiently  utilize  the  available  spectrum.  The  unprecedented 
radio  agility  envisioned,  calls  for  fast  and  accurate  spectrum 
sensing  over  a  wide  bandwidth,  which  challenges  traditional 
spectral  estimation  methods  typically  operating  at  or  above 
Nyquist  rates.  Capitalizing  on  the  sparseness  of  the  signal 
spectrum  in  open-access  networks,  this  paper  develops  com¬ 
pressed  sensing  techniques  tailored  for  the  coarse  sensing  task 
of  spectrum  hole  identification.  Sub-Nyquist  rate  samples  are 
utilized  to  detect  and  classify  frequency  bands  via  a  wavelet- 
based  edge  detector.  Because  spectrum  location  estimation 
takes  priority  over  fine-scale  signal  reconstruction,  the  pro¬ 
posed  novel  sensing  algorithms  are  robust  to  noise  and  can 
afford  reduced  sampling  rates. 

Index  Terms —  spectrum  estimation,  compressed  sens¬ 
ing,  sub-Nyquist  sampling,  wavelet  transform,  cognitive  radio 

1.  INTRODUCTION 

The  emerging  paradigm  of  Dynamic  Spectrum  Access  shows 
promise  to  alleviate  today’s  spectrum  scarcity  problem  by  ush¬ 
ering  in  new  forms  of  spectrum  agile  wireless  networks  [1]. 
Key  to  this  new  paradigm  are  cognitive  radios  (CRs)  that  are 
aware  of  and  can  sense  the  environments,  and  perform  func¬ 
tions  to  best  serve  their  users  without  causing  harmful  interfer¬ 
ence  to  other  authorized  users  [2].  As  such,  the  first  cognitive 
task  preceding  any  form  of  dynamic  spectrum  management 
is  to  develop  wireless  spectral  detection  and  estimation  tech¬ 
niques  for  sensing  and  identification  of  available  spectrum. 

Spectrum  sensing  in  the  wideband  regime  faces  consider¬ 
able  technical  challenges.  The  radio  front-end  can  employ  a 
bank  of  tunable  narrowband  bandpass  filters  to  search  one  nar¬ 
row  frequency  band  at  a  time.  In  each  narrowband,  existing 
spectrum  sensing  techniques  perform  either  energy  detection 
[3]  or  feature  detection  [2].  It  requires  an  unfavorably  large 
number  of  RF  components  and  the  tuning  range  of  each  fil¬ 


ter  is  preset.  Alternatively,  a  wideband  circuit  utilizes  a  sin¬ 
gle  RF  chain  followed  by  high-speed  DSP  to  flexibly  search 
over  multiple  frequency  bands  concurrently  [4].  A  major  im¬ 
plementation  challenge  lies  in  the  very  high  sampling  rates 
required  by  conventional  spectral  estimation  methods  which 
have  to  operate  at  or  above  the  Nyquist  rate.  Meanwhile,  due 
to  the  timing  requirements  for  rapid  sensing,  only  a  limited 
number  of  measurements  can  be  acquired  from  the  received 
signal,  which  may  not  provide  sufficient  statistic  when  tradi¬ 
tional  linear  signal  reconstruction  methods  are  employed. 

This  paper  aims  at  fast  spectrum  sensing  at  affordable  com¬ 
plexity.  A  couple  of  key  premises  are  capitalized  to  alleviate 
the  stringent  sampling  requirements  in  the  wideband  regime. 
First,  we  take  a  multi-resolution  approach  to  decompose  the 
cognitive  sensing  task  into  two  stages.  The  first  stage  is  coarse 
sensing  to  detect  non-overlapping  spectrum  bands  and  classify 
them  into  black,  gray  or  white  spaces,  depending  on  whether 
the  power  spectral  density  (PSD)  levels  are  high,  medium  or 
low  [2].  Based  on  the  spectrum  sharing  mechanism  adopted 
[1],  the  second  stage  of  fine-scale  spectral  shape  estimation 
is  performed  only  when  needed,  and  mostly  confined  within 
the  available  (narrowband)  white  spaces  to  alleviate  the  sam¬ 
pling  requirements.  Second,  we  recognize  that  the  wireless 
signals  in  open-spectrum  networks  are  typically  sparse  in  the 
frequency  domain.  This  is  due  to  the  low  percentage  of  spec¬ 
trum  occupancy  by  active  radios  -  a  fact  motivating  dynamic 
spectrum  management.  For  sparse  signals,  recent  advances  in 
compressed  sensing  have  demonstrated  the  principle  of  sub- 
Nyquist-rate  sampling  and  reliable  signal  recovery  via  com¬ 
putationally  feasible  algorithms  [5,  6,  7,  8]. 

Tailored  to  the  above  distinct  nature  of  CR  sensing,  this 
paper  derives  novel  compressed  sensing  algorithms  for  the 
coarse  sensing  task  of  spectrum  band  classification.  Random 
sub-Nyquist-rate  samples  are  employed  to  formulate  an  op¬ 
timal  signal  reconstruction  problem,  which  incorporates  the 
wavelet-based  edge  detector  we  recently  developed  in  [10]  to 
recover  the  locations  of  frequency  bands.  Because  spectrum 
location  estimation  takes  priority  over  fine-scale  signal  recon¬ 
struction,  the  novel  sensing  algorithms  are  robust  to  noise  and 
can  afford  reduced  sampling  rates. 
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2.  SIGNAL  MODEL  AND  PROBLEM  STATEMENT 

Suppose  that  a  total  of  B  Hz  in  the  frequency  range  [/o,  /at]  is 
available  for  a  wideband  wireless  network.  A  CR  receives  the 
signal  r[t)  that  occupies  N  consecutive  spectrum  bands,  with 
their  frequency  boundaries  located  at  /o  <  /i  <  •  •  •  /tv-  The 
frequency  response  of  r[t)  is  illustrated  in  Fig.  1.  Depending 
on  whether  the  PSD  level  is  high,  medium  or  low,  each  fre¬ 
quency  segment  can  be  categorized  into  black,  gray  or  white 
spectrum  spaces  [2].  White  holes,  and  sometimes  gray  spaces, 
can  be  picked  by  the  CR  for  opportunistic  transmission,  while 
the  black  holes  are  to  be  avoided  for  interference  control. 


Fig.  1.  N  frequency  bands  with  piecewise  smooth  PSD. 

Suppose  that  the  time  window  for  sensing  is  f  G  [0,  MTq], 
where  Tq  is  the  Nyquist  sampling  rate.  Using  Nyquist  sam¬ 
pling  theory,  M  samples  are  needed  to  recover  r[t)  without 
aliasing.  A  digital  receiver  converts  the  continuous-domain 
signal  r{t)  to  a  discrete  sequence  x*  G  of  length  K.  The 
sampling  process  can  be  expressed  in  discrete-time  domain  in 
the  following  general  form; 

xt  =  S^r*  (1) 

where  S  is  an  M  x  IT  projection  matrix  and  r*  represents 
the  M  X  1  vector  with  elements  rt[n]  =  r{t)\t=nTo,  ^  = 

1, . . . ,  M.  Columns  of  S  can  be  viewed  as  a  set 

of  basis  signals  or  matched  filters,  while  the  measurements 
{^tW\}k=i  essence  the  projection  of  r{t)  onto  the  ba¬ 

sis.  The  model  in  (1)  subsumes  all  sampling  schemes  yield¬ 
ing  linear  measurements.  For  example,  S  =  Itw  represents 
Nyquist-rate  uniform  sampling,  where  I  tu  is  the  size-M  iden¬ 
tity  matrix;  S  =  Fm  amounts  to  frequency-domain  sampling, 
where  Fm  is  the  M-point  unitary  discrete  Fourier  transform 
(FT)  matrix.  When  K  <  M,  reduced-rate  sampling  arises. 
We  focus  on  non-adaptive  measurements  where  S  is  preset. 

The  goal  of  CR  sensing  is  to  classify  and  estimate  the 
spectrum  of  r[t)  given  the  sample  set  xj,  where  K  <  M 
is  possible.  Spectrum  classification  refers  to  identifying  the 
number  of  subbands  N  and  their  locations  {[/,,  /i-i-i]}^o^, 
and  classifying  them  into  black,  gray  or  white  spaces.  Spec¬ 
trum  estimation,  on  the  other  hand,  can  have  different  objec¬ 
tives:  either  to  estimate  the  frequency  response  of  r{t)  within 
the  entire  wideband,  or  confine  the  estimation  to  be  within  the 
identified  (narrowband)  white  spaces  only.  This  paper  primar¬ 
ily  concerns  the  coarse  sensing  task  of  spectrum  classification. 


3.  MULTI-STEP  COMPRESSED  SENSING 

Our  first  approach  to  reduced-complexity  spectrum  sensing 
takes  the  following  four  steps:  i)  compressed  random  sam¬ 
pling  to  generate  measurements  xj  from  r(f);  ii)  reconstruc¬ 
tion  of  the  frequency  response  y  j  =  FmYi  from  x*;  Hi)  esti¬ 
mation  of  frequency  band  number  N  and  locations 
based  on  r  / ;  and,  iv)  estimation  of  the  average  amplitude  of 
r  f  within  each  identified  band  for  spectrum  classification.  It 
is  worth  emphasizing  that  Step  ii)  recovers  the  accurate  fine- 
resolution  signal  spectrum  r/  represented  by  M  frequency 
samples  at  the  Nyquist  rate,  while  the  available  measurement 
set  Xt  is  of  a  reduced  size  of  K{<  M)  elements. 

3.1.  Sub-Nyquist-rate  Sampling 

Let  IF  denote  the  non-zero  frequency-domain  support  of  r{t) 
in  the  noise-free  case.  In  open-spectrum  networks,  it  gener¬ 
ally  holds  that  \F\  B  [4],  indicating  the  sparseness  nature 
of  r{t).  Equivalently  speaking,  the  M  x  1  frequency  response 
vector  r/  contains  on  average  K},  :=  non-zero 

elements  when  noise  free,  and  Kf,  M.  The  key  results 
in  compressed  sensing  stated  that  the  sparse  vector  r  j  can  be 
recovered  asymptotically  from  K  (<  M)  samples  of  r{t),  as 
long  as  K  >  Kb.  These  samples  xt  can  be  generated  from  (1) 
via  universal  non-uniform  sampling  [5]  or  random  sampling 
[6,  7],  both  of  which  can  enable  perfect  recovery  of  r  /  when 
free  of  noise.  To  distinct,  we  denote  a  reduced-rate  sampling 
matrix  as  Sc  of  dimension  M  x  K,  where  Kb<K  <^M.  A 
simple  example  of  Sc  is  a  selection  matrix  that  randomly  re¬ 
tains  K  columns  of  the  size-M  identity  matrix,  which  means 
that  K  —  M  time  instants  on  the  sampling  grid  are  skipped. 

3.2.  Spectrum  Reconstruction 

With  the  K  measurements  xj  =  S^r*,  we  now  estimate  the 
frequency  response  of  r{t)  in  the  form  of  r  y  =  Fj^r.  For 
a  given  linear  sampler  Sc  :  ,  we  seek  a  nonlinear 

reconstruction  function  R{-)  :  that  offers  an  ap¬ 
proximate  reconstruction  of  r  /  G  from  xj  G  based  on 

the  linear  transformation  equality  Xt  =  {S'^F^  )r/;c.f.  (1). 
This  is  a  linear  inverse  problem  with  sparseness  constraint, 
which  is  NP-hard.  A  conceptually  intuitive  approach  to  signal 
reconstruction  is  the  Basis  Pursuit  (BP)  technique  [9],  which 
transforms  the  sparseness  constraint  on  r  /  into  a  convex  opti¬ 
mization  problem  solvable  by  linear  programming: 

f/ =  arg  min||r/||i,  s.t.  (S^Fjj/jr/ =  x*.  (2) 

!■/ 

Besides  BP,  a  number  of  efficient  reconstruction  methods  ex¬ 
ist,  including  orthogonal  matching  pursuit  (OMP)  algorithm 
and  tree-based  OMP  (TOMP)  algorithm  [8].  Since  the  mea¬ 
surements  can  be  complex- valued,  we  find  it  convenient  to  use 
TOMP  in  our  simulations,  but  for  illustration,  formulate  our 
signal  reconstruction  problem  based  on  BP,  as  in  (2). 
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3.3.  Band  Location  Estimation 

Having  estimated  r/,  we  turn  to  the  wavelet-based  edge  de¬ 
tector  in  [10]  for  detecting  the  number  and  frequency  loca¬ 
tions  of  spectrum  spaces.  The  basic  idea  is  to  view  the  entire 
wideband  under  scrutiny  as  a  train  of  consecutive  frequency 
subbands,  where  the  PSD  is  smooth  within  each  subband,  but 
exhibits  a  discontinuous  change  between  adjacent  subbands. 
These  irregularities  are  in  fact  edges  in  PSD,  which  carry  key 
information  on  the  locations  and  intensities  of  spectrum  holes. 
To  further  simplify  and  expedite  the  coarse  sensing  stage,  we 
approximately  treat  the  spectral  amplitudes  within  each  sub¬ 
band  to  be  almost  flat,  at  an  unknown  level  a  „  over  the  n-th 
band.  These  modeling  approximations  are  invoked  to  reduce 
the  overall  wideband  sensing  complexity.  If  needed,  the  sens¬ 
ing  quality  can  be  refined  after  spectrum  holes  are  identified. 

Based  on  these  modeling  assumptions  and  with  reference 
to  Fig.  1,  wideband  sensing  can  be  viewed  as  an  edge  detec¬ 
tion  problem  in  an  image  depicted  by  f  /  in  frequency.  Edges 
in  this  image  correspond  to  the  locations  of  frequency  discon¬ 
tinuities  which  are  to  be  identified.  The  wavelet 

approach  is  well  motivated  for  edge  detection  [11].  Given  ry, 
we  re-cast  the  edge  detector  in  [10]  in  discrete  form. 

Let  (f>{f  )  be  a  wavelet  smoothing  function  with  a  compact 
support.  The  dilation  of  0(/)  by  a  scale  factor  s  is  given  by 

For  dyadic  scales,  s  takes  values  from  powers  of  2,  i.e.,  s  = 
V,j  =  1,2,...,J.  Let$,(r)  :=  J-i{0«(/)}  =  $(sr) 
represent  the  inverse  FT  of  the  wavelet  function.  The  contin¬ 
uous  wavelet  transform  of  R{f  ){^  r y)  is  given  by  [10] 

WMf)  =  ^{>V.r(T)}  =  ■  $(st)}.  (4) 

Mapping  WsR{f  ),  r(T)  and  <1>(st)  to  their  length-M  discrete 
counterparts  y^,  r*  and  <l>s  respectively,  (4)  is  equivalent  to 

ys  =  FMdiag{$s}rt.  (5) 

Replacing  rt  in  (5)  by  its  estimate  rt  =  Fjj^^ry,  we  reach 
the  estimated  wavelet  transform 

y*  =  FMdiag{<I>jFj[:/ry.  (6) 

The  derivative  wavelet  of  r  y  at  scale  s  is  given  by  z,,  with 
elements  {zs[n]}^^i  in  the  form 

s-^  {WsR{f))  ^  Zs  :  Zs[n\  =  -  ys[n  -  1].  (7) 

The  boundaries  {fn}n=i  thus  be  acquired  by  picking  the 
local  maxima  of  the  wavelet  modulus  z  g ,  while  the  band  num¬ 
ber  N  is  determined  by  the  number  of  local  peaks  [10,  11]. 


3.4.  Frequency  Response  Amplitude  Estimation 

The  estimated  boundaries  {fn}n=i  correspond  to  A  —  1  se¬ 
lected  indices  {/„  :  /„  =  /o  -I-  /„A,A  =  B/M}  in  the 
frequency  response  vector  r  y .  Elements  of  r  y  between  a  pair 
of  adjacent  indices  belong  to  the  same  frequency  band.  The 
average  frequency  response  amplitude  of  the  n-th  band, 
can  thus  be  computed  as 

1 

- TT  n=l,...,N.  (8) 

Jn— 1  “r  -L  .  , 

1  =  1^-! 

This  simple  and  coarse  estimator  in  (8)  allows  us  to  categorize 
the  detected  frequency  bands  into  black,  gray,  or  white  spaces 
[2],  depending  on  whether  {ctn}  are  high,  medium  or  low. 

4.  ONE-STEP  COMPRESSED  SENSING 

To  further  reduce  the  implementation  complexity  of  coarse 
spectrum  sensing,  we  now  ask:  can  we  directly  detect  and 
estimate  the  frequency  band  locations  from  the  compressed 
measurements  xt  in  (1),  without  having  to  recover  the  detailed 
frequency  response  ry  ?  We  address  this  question  by  deriving 
signal  recovery  formulation  for  wavelet-based  edge  detection. 

Recall  from  Section  3.3  that  the  band  locations  can  be  re¬ 
covered  from  the  {N—1)  peaks  of  the  derivative  wavelet  mod¬ 
ulus  Zg  e  .  When  N  ^  M  (which  is  generally  the  case), 
Zg  can  be  treated  as  a  sparse  vector,  with  only  a  few  non¬ 
trivial  elements  located  at  frequency  band  boundaries;  c.f.,  see 
Fig.  3  for  graphical  validation.  As  such,  z  g  can  be  recovered 
under  the  sparseness  constraint,  provided  that  we  can  find  a 
linear  transformation  equality  linking  Zg  to  the  compressed 
measurement  vector  x*. 

To  this  end,  we  rewrite  (7)  in  matrix- vector  form  as  z  g  = 
r Ys ,  where  F  is  the  differentiation  matrix  given  by 

r  1  0  •••  0 1 


0 

-  ®  ^  -  MxM 

Putting  (9)  and  (5)  together,  we  obtain: 

D  =  (Fm  diag{$g})  ^  yg  =  (Fj;/diag{<I>g})  ^  •  F'^  -Zg. 

' - V - " 

:=G 

(10) 

Noting  that  x*  =  S^Tr*,  and  that  Zg  is  sparse,  we  reach  the 
following  BP-based  optimization  formulation: 

Zg  =  arg  min  ||zg||i,  s.t.  xt=(SjG)zg.  (11) 

Zs 

Subsequently,  the  band  boundaries  {/„}  can  be  acquired 
from  the  locations  of  those  non-zero  elements  in  z  g ,  obviating 
the  involved  step  of  frequency  response  estimation  on  r  y. 
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5.  SIMULATIONS 

We  consider  a  wide  band  of  interest  in  the  range  of  /g  + 
[50, 150]A  Hz,  where  A  is  the  frequency  resolution.  Fig.  2 
illustrates  the  spectral  amplitude  |i?(/)|  observed  by  a  CR. 
During  the  observed  burst  of  transmissions  in  the  network, 
there  are  a  total  of  A  =  6  bands,  with  frequency  boundaries 
at  {/n}Li  =  /o  +  [60, 68, 83, 119, 123, 150]A  Hz.  Among 
these  bands  (marked  in  Fig.  2),  Bi,  and  B^  have  rela¬ 
tively  high  signal  amplitude  at  levels  16,  20,  and  24,  respec¬ 
tively,  while  B2  has  low  amplitude  at  a  level  of  2.  The  rest  two 
bands,  B4  and  Be  are  not  occupied  and  are  thus  white  spec¬ 
trum  holes.  The  sampling  lower  bound  is  thus  K t /M  «  40%. 

For  compressed  sensing,  the  compression  ratio  K/AI  is 
set  to  vary  from  50%  to  100%  with  reference  to  the  Nyquist 
rate.  The  noise  level  is  =  8  dB.  The  sampler  S  =  Sc  used 
in  (1)  is  uniformly  random.  Fig.  2  indicates  that  the  signal 
recovery  quality  (via  TOMP)  improves  as  K/M  increases. 

In  the  wavelet-based  edge  detector,  Gaussian  wavelets  are 
used  at  four  dyadic  scales  s  =  2^ ,  j  =  1,2,3, 4.  Fig.  3  de¬ 
picts  the  multiscale  wavelet  products  computed  from  (7)  [10]. 
Edges  in  the  R{f)  are  clearly  captured  hy  the  wavelet  trans¬ 
form  in  all  curves.  As  the  scale  factor  increases,  the  wavelet 
transform  becomes  smoother  within  each  frequency  band,  re¬ 
taining  the  lower- variation  contour  of  the  noisy  PSD. 

For  frequency  hand  location  estimation.  Fig.  4  depicts 
the  normalized  root  mean-square  estimation  errors  (RMSE) 

B~^ \  fn  —  /nP  with  respect  to  both  the  compres¬ 
sion  ratio  K/M  and  the  inverse  noise  level  n“^.  When  either 
the  number  of  samples  is  very  small  or  the  noise  is  very  strong, 
there  exhibits  an  estimation  error  floor.  Nevertheless,  the  at¬ 
tained  degree  of  estimation  accuracy  is  beneficial  to  effecting 
CR  agility  at  affordable  sampling  cost.  Robustness  to  sample 
quantization  errors  is  also  illustrated. 
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Fig.  2.  signal  frequency  response:  (top)  noise-free  [Ay  |;  (rest) 
recovered  |  A/|  at  compressing  ratios  K/M  =  50%,  75%,  1. 


Fig.  3.  wavelet  modulus  for  edge  detection  at  scales  s  =  2P 


Fig.  4.  frequency  band  location  estimation  errors  RMSE. 
Strong  compression  of  K/M  =  33%  is  demonstrated.  Im¬ 
pact  of  quantization  (bits/sample)  is  shown  for  K/M  =  50%. 
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